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Determination  of  the  Aurorai  Ovai  Q 
Index  from  the  Air  Weather  Service  K 

index 


1.  INTRODncnON 

Ionospheric  models  such  as  the  ITS-78  (Barghausen  et  al^  1969)  and/or  lONCAP  (Lloyd  et 
al,^  1978)  are  unable  to  predict  the  ionospheric  parameters  accurate^:  the  critical  frequency,  the 
height  of  the  peak  electron  density,  and  the  thickness  of  the  layers,  for  the  high  latitude  region 
affected  by  auroral  enhancements,  and  for  the  midlatitude  F  layer  trough,  vdiich  is  the  region  of 
ionospheric  depletion.  To  overcome  this  deficiency,  the  Air  Force  Global  Weather  Central 
(AFGWC)  developed  the  Polar  Ionospheric  model  (see  Millman  et  al®)  by  modifying  the  lTS-78 
model  for  the  auroral  enhancement  and  for  the  F-layer  trough,  to  support  qrstems  like  Over  the 
Horizon  Backscatter  (OTHB)  radars  operating  in  the  mid  and  high  latitude  regions. 

The  input  for  the  rrS78/IONCAP  models  are  the  Julian  day.  the  universal  time,  and  the 
simspot  number.  The  AWS  uses  the  effective  sunspot  number  SSNeff  determined  firom  the  f^2 
observations  firom  the  AWS  network  of  the  Digital  Ionospheric  Sounding  Systems  (DISS)  and  other 
ionosondes.  In  addition  to  the  SSNeff.  the  Polar  Ionospheric  Model  needs  the  magnetic  activity 
indices:  the  planetary  index  Kp  and  the  auroral  oval  index  Q.  which  control  respectively  the 
midlatitude  trough  and  the  auroral  oval  algorithms  in  the  model.  For  generating  the  real  time  'Kp' 
index  (AWS-IQ  the  AWS  maintains  a  network  of  eight  magnetometer  stations  covering  the 
European  and  the  American  sectors  (see  Prochaska^  auid  Dandekar®).  At  present  the  auroral  Q 
index  is  measured  at  a  single  station  at  Sodankyla.  Finland.  This  Q  index  is  not  available  in  real 
time  to  the  AWS. 

Ideally  the  Q  index  should  be  measured  at  an  auroral  latitude  on  its  magnetic  midnight 
meridian.  When  the  Q  index  is  measured  at  a  single  station,  the  observed  value  is  lower*  when 
the  rotation  of  the  earth  moves  the  station  away  firom  the  auroral  oval.  The  AWS  has  no  network 
similar  to  that  for  the  'Kp'  (AWS-K)  index,  for  measuring  the  Q  index.  Because  the  Q  index  is  not 
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available  in  real  time,  Hardy^  (private  communication)  developed  a  scheme  of  using  the  DMSP- 
satellite  based  particle  precipitation  measurements.  On  the  basis  of  the  analysis  of  the  DMSP 
data^-^  (Gussenhoven  et  al  1981,  Gussenhoven  et  al  1983),  the  satellite  measurements  provide  the 
instantaneous  Kp  index.  Hardy^  (private  communication)  derived  an  algorithm  to  relate  the  Kp 
and  Q  indices  by  matching  the  Starkov  Q=0  optical  auroral  oval*°  with  the  Kp=0  particle 
precipitation  oval.  During  the  OTH-B  tests  we  noticed  that  the  Q  index  determined  from  the 
Hardy  algorithm  is  high,  especially  for  low  Kp  values.  This  report  gives  a  correction  for  the 
problem. 

The  first  section  discusses  the  Hardy  algorithm.  The  second  section  seeks  other  sources  for 
determining  a  reliable  relation  between  the  Q  and  Kp  indices.  A  relation  more  precise  than  the 
Hardy  algorithm  will  provide  improved  predictions  for  the  aiunral  oval  phenomena.  The  auroral 
oval  is  one  of  the  important  clutter  sources  for  the  systems  like  the  OTHB-ECRS  and  WCRS 
radars,  for  which  their  northern  sector  coverage  of  high  latitude  region  is  affected  by  the 
ionospheric  irregularities  associated  with  the  auroral  oval  and  the  trough.  The  improved 
prediction  of  the  auroral  oval  will  help  in  relating  the  high  latitude  clutter  more  closely  with  the 
auroral  oval  morphology. 

2.  BACKGROUND 

In  the  absence  of  the  real  time  auroral  oval  Q  index,  the  AWS  uses  real  time  particle 
precipitation  boundaries  observed  by  satellite  (Defense  Meteorological  Satellite  Program  DMSP) 
borne  sensors.  These  polar  orbiting  satellites  provide  a  capability  for  determining  the  precipitation 
boundary  every  1 10  minutes.  Gussenhoven  et  al®-®  derived  hourly  empirical  relations  between  the 
precipitation  boundary  and  the  planetary  magnetic  activity  index  Kp.  By  using  these  relations  Kp 
can  be  determined  from  the  DMSP  particle  precipitation  data.  For  convenience  the  latest  set  of 
the  coefiicients  for  the  straight  line  fits  available  from  the  Madden  and  Gussenhoven  study  is 
reproduced  in  Table  1  (Table  3  from  their  report).  The  table  presents  the  results  of  the  least 
squared  deviation  straight  line  fits  between  the  observed  boimdary  ’Lp’  and  the  corresponding  Kp 
index  for  hourly  intervals  of  the  corrected  magnetic  local  time  (CGT).  The  first  column  gives  the 
time  interval  T  (in  CGT)  and  the  second  column  lists  the  number  of  observations.  The  next  three 
columns  provide  the  slope  Ap,  the  intercept  Bj.  (at  Kp=0),  and  the  correlation  coefficient  'C  for  the 
fit.  The  correlation  coefficient  provides  the  level  of  confidence  for  the  empirical  fit.  The  percent 
population  closefy  following  the  fit  is  between  and  C.  Thus  with  ICI=0.707,  50  percent  to  71 
percent  (more  than  half)  of  the  population  of  the  data  follows  the  fit.  A  look  at  the  last  column  of 
Table  1  shows  a  lack  of  correlation  ( I C I  <0,6)  for  the  intervals  of  2-4  and  of  9-16  hours.  When  the 
correlation  is  significant,  the  observed  boundary  can  be  used  for  estimating  the  Kp  index  by  using 
the  equation. 


Lp  =  Ap  Kp  +  Bp 


(1) 


2 


Table  1.  Regression  Coe£Bcients  of  Auroral  Equatorward  Boundaries  vs 
(Madden-Gussenhoven  Data  1990) 


REGRESSION  COEFTICIENTS  OF  AURORAL 
EQUATORWARD  BOUNDARIES  VS  KP 
ALL  BOUNDARIES,  1983-1990 


Magnatic 
Local  Tine 

Miater  of. 
Boundaries 

Slope 

Bt 

00-01 

2349 

67.57 

-1.62 

01-02 

43 

68.31 

-1.39 

02-03 

41 

68.92 

03-04 

2741 

67.07 

-1.50 

04-05 

12900 

66.56 

-1.82 

05-06 

20682 

67.28 

-1.79 

06-07 

16186 

68.10 

-1.83 

07-08 

13369 

68.17 

-1.65 

08-09 

17422 

68.98 

-1.58 

09-10 

17873 

69.13 

-1.24 

10-11 

6565 

68.99 

-1.00 

11-12 

2218 

68.54 

12-13 

1056 

68.90 

13-14 

1296 

70.76 

14-15 

1564 

71.48 

■Hi 

15-16 

1785 

72.73 

-1.12 

16-17 

3332 

73.22 

-1.46 

17-18 

8193 

72.20 

-1.48 

18-19 

19946 

71.64 

-1.64 

19-20 

17347 

71.09 

-1.85 

20-21 

17539 

69.71 

-1.66 

21-22 

14175 

69.25 

-2.07 

22-23 

16234 

67.89 

-1.88 

23-24 

11205 

67.18 

-1.75 

where  is  the  observed  boundary  in  CG  latitude  at  0ven  CGT  T, 

Af  is  the  slope  of  the  fit  for  time  T. 

is  the  intercept  of  the  fit  for  time  T.  and 
Kp  is  the  level  of  magnetic  activity  at  time  T. 

The  intercept  Bj.  is  related  to  the  radius  of  the  auroral  oval  for  KpsQ,  and  the  slope  is  the 
increase  in  the  radius  of  the  auroral  oval  for  each  unit  increase  in  I4>- 

Feldstein  and  Starkov^^  determined  the  auroral  oval  boundary  from  All  Slqr  Camera  (ASCA) 
observations,  and  related  it  to  the  oval  index  Q  measured  on  the  magnetic  midnight  meridian.  We 
will  refer  to  these  as  Starkov  optical  auroral  ovals.  In  the  nightside  the  particle  precipitation 
boundary  is  practically  coincident  with  the  optical  auroral  oval  boundary.  Therefore  Heudy^ 
matched  the  midnight  boundary  of  the  DMSP  (Kp=0)  particle  precipitation  with  the  Starkov  (Q=0) 
optical  auroral  ovals  for  determining  the  Kp-Q  relation.  The  Starkov  optical  auroral  oval 
boundary*®  is  given  by 


At  =  90-1(18  +  0.9Q  +  5. 1  cos  (15T-12))  for  Q^tO. 


(2) 


The  resulting  empirical  relation  derived  by  Hardy  is  given  by 


9  =  2KP-0.35 


(3) 


In  his  report  Hardy  (private  communication)  notes  that  the  midnight  precipitation  and  optical 
auroral  boundaries  match,  but  the  difference  in  the  respective  boundaries  is  successively  larger 
as  one  moves  away  fiom  the  midnight  sector  to  the  noon  sector.  He  reports  that  at  Kp=0  the 
noon  latitude  of  the  Starkov  optical  auroral  oval  is  7.8®  poleward  from  the  particle  precipitation 
oval,  and  for  1^=3  this  latitude  difference  is  5®. 

Figure  1  shows  the  plot  of  Kp-Q  data  obtained  for  52  days  (381  observations)  during  the  tests 
of  the  ECRS  and  the  WCRS  radars.  About  76  percent  of  the  observations  cluster  around  the  line 
given  by  Eq.  (3).  Our  inquiry  to  the  AWS  showed  that  when  satellite  data  are  absent,  the  ground 
based  Kp(AWS-K)  data  are  put  into  the  Hardy  Ek}.  (3)  to  determine  the  Q  indeic.  In  Figure  1  these 
ground  based  discrete  data  levels  lie  right  on  the  straight  line  given  by  Ek].  (3)  and  produce  the 
clusters.  In  this  small  data  set,  the  nonclustered  data  show  that  less  than  24  percent  of  the 
observations  are  satellite  based.  We  feel  that  in  Figure  1,  the  computed  Q  index  is  overestimated, 
especially  for  Kpsl  and  2  (Qsl.65  and  3.65  respectivety),  because  experience  shows  that  the  9=3 
condition  is  approximately  equal  to  that  of  i4>=3. 

Not  only  is  the  point  Kp=3,  Q=3  below  zmy  of  the  lines;  it  is  well  below  the  average  of  the 
satellite  data.  However,  the  time  of  the  satellite  reading  (in  CGI)  is  not  available.  Without  this 
time  information,  the  movement  of  the  auroral  precipitation  oval  with  respect  to  the  optical  oval 
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Determination  of  the  Q  Index  Using  Hardy  Algorithm  (Equation  3) 


must  remain  unknown.  In  addition,  the  sample  of  satellite  data  is  not  large  enou^  to  be  reliable 
statistically.  These  limitations  motivated  a  search  for  a  relationship  that  would  reflect  what  we  do 
know  about  the  1^0  relationship,  including  the  fact  that  Q=3  for  Kps3. 

In  the  following  we  will  compare  the  Starkov  Q  (optical)  and  the  Kp  particle  precipitation 
ovals,  to  study  the  S3rstematic  differences  between  the  Q  and  Kp  auroral  ovals  for  the 
determination  of  a  more  appropriate  Kp-Q  relation.  Hardy  used  the  hourly  equations  from  Table  1 
for  detennining  the  Kp  index.  In  this  study,  two  least  squared  deviation  empirical  circles  were 
fitted,  one  to  A,,  and  the  other  to  Bj..  For  this  purpose  onty  the  Af-Bj  values  with  I C I  =>0.6  from 
Table  1  were  used.  TTie  equation  derived  for  the  colatitude  of  the  intercept  is 

B't=  90-Bt.=  20.27  +  3.37  cos(15T-41.8)  (4) 

and  the  equation  for  the  slope  is 

Ar=  1.69  +  0.24  cos(15T-13.8)  (5) 

In  each  of  the  above  equations  the  first  term  on  the  RHS  gives  the  radius  of  the  circle,  the 
constant  multiplier  of  the  cosine  term  provides  the  offset  (in  latitude  from  the  CG  pole)  of  the 
circle  from  the  center,  and  the  number  in  the  bracket  of  the  cosine  term  presents  the  angular 
deviation  (hour*  15)  from  the  midnight  meridian.  Combining  the  two  equations  provides  an 
empirical  fit  between  the  observed  boundary  and  the  Kp  index  for  the  particle  precipitation  ovals. 
Substituting  E^s.  (4)  and  (5)  in  E^.  (1)  we  get 

At  =  90  -  (20.27  +  3.37cos(15T-42)  +  Kp  {-[1.69  +  0.24cos(15T  -  14)]})  (6) 

This  single  equation  (6)  can  be  used  in  place  of  16  T  (time)  sets  of  Eq.  (1)  from  Table  1. 

The  comparison  of  Ekjs.  (2)  and  (6)  is  summarized  in  Table  2.  The  table  shows  that  the  radius 
of  the  Kp=0  particle  precipitation  oval  is  0.6°  larger  than  that  of  the  9=0  optical  auroral  oval.  The 
center  of  the  Q  optical  auroral  oval  is  1.7°  equatorward  of  that  of  the  Kp  auroral  oval.  The  radius 
of  the  Kp  auroral  oval  increases  at  1.7°  per  unit  Kp.  whereas  the  radius  of  the  optical  auroral  oval 
increases  at  0.9°  per  unit  Q.  Thus,  the  respective  a-iroral  ovals  have  different  starting  locations 
and  different  rates  of  expansion  with  the  increase  in  the  level  of  magnetic  activity.  Therefore 
forcing  a  match  between  the  two  auroral  ovals  at  the  midnight  meridian  results  in  increasing 
differences,  with  increasing  level  of  magnetic  activity,  in  their  locations  at  meridians  away  from 
the  (CG)  midnight. 

The  data  from  Table  1  and  the  empirical  circle  fit  from  Ekj.  (6)  are  shown  in  Figure  2.  The  data 
points  with  I C I  <0.6  eu'e  not  used  in  the  study  and  are  not  plotted  in  Figure  2.  The  cause  of  the 
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Figure  2.  Least  Squared  Deviation  Circle  Ht  to  Madden-Gussenhoven  Data  (Equation  6) 


systematic  deviation  of  00-01-02  hours'  data  firom  the  experimental  fit  is  not  known. 

In  the  following  we  compare  the  Q  optical  auroral  ovals  with  the  Kp  pzutlcle  precipitation 
ovals  (using  Eq.  (6))  for  checking  the  validity  of  the  Kp-Q  relation  of  Eq.  (3). 

In  addition  to  Eq.  (2),  Starkov  provides  a  separate  equation*®  for  Q=0  (For  the  polar  ionospheric 
model  the  AFGWC  uses  the  Feldstein-Starkov  (1967)  oval  **  for  Q=0  and  the  Starkov  equation  (Ekj. 
(2)  shown  above)  for  Q^).  Hie  equation  is 


hj.=  0.  -  [(72.5-3.4  cos(l5T-18)+0.7cos(2*T*15-45)]. 


(7) 


Starkov  optical  auroral  ovals,  for  Q=0.  fipom  Eqs.  (2)  and  (7)  are  shown  in  Figure  3.  Note  that  the 
two  Q=0  auroral  ovals  have  the  same  boundary  at  16-18  CGT  and  06-07  CGT,  and  the  difference 
is  largest  (2.9°)  around  00-02  CGT,  with  the  Eq.  (2)  auroral  oval  more  equatorward  than  that  finm 
Eq.  (7).  and  for  09-10  CGT,  the  auroral  oval  from  Eiq.  (7)  is  1.2°  equatorward  of  that  frnm  Eq.  (2). 
For  matching  the  midnight  boimdary  of  the  auroral  oval  of  Ek).  (2)  with  that  from  Eq.  (7)  we  need 
Q=-3.  To  match  noon  boundaries  we  need  Q=1  for  Eq.  (2).  Thus,  there  is  a  significant  difference 
between  the  Q=0  auroral  ovals  from  Ekjs.  (2)  and  (7). 

The  next  figure  (Figure  4)  compares  auroral  ovals,  optical  for  Q=0  (Eq.  (3))  and  particle 
precipitation  for  Kp=0  (Eq.  (7)).  Note  that  both  the  auroral  ovals  match  for  20-01  CGT  but  for  10- 
12  CGT  the  DMSP  Kp  particle  precipitation  auroral  oval  is  4.5°  equatorward  of  the  optical  auroral 
oval. 

Figure  5  shows  the  optical  auroral  oval  for  9=9,  with  particle  precipitation  ovals  for  Kp=2  and 
4.  Note  that  for  Kp=4  and  9=9  the  auroral  oval  boundaries  match  at  01-02  CGT,  whereas  at  10- 
1 1  CGT  the  Q=9  boundary  is  closer  to  that  for  Kp=2. 

Thus,  we  see  that  there  is  no  good  correlation  between  the  Starkov  optical  and  the  particle 
precipitation  auroral  ov2ds,  especially  in  the  day  sector.  Significant  differences  exist  between  the 
optical  9  and  the  particle  precipitation  Kp  ovals,  because  the  DMSP  boundaries  are  from  the 
particle  precipitation  data  with  a  lower  threshold  of  500  eV  that  is.  with  higher  sensitivity,  than 
those  for  the  optical  emission  boundaries  for  the  Feldstein-Starkov  auroral  ovals  determined  fi'om 
the  All  Slqr  Camera  Systems  (ASCA),  which  had  lower  sensitivity  because  of  its  higher  threshold  of 
1-2  keV.  Whalen  et  al.  (1971)  observed  6300  A  soft  particle  precipitation  equatorwards  of  the 
auroral  E  layer  in  the  noon  sector.*®  Winningham  and  Heikkila  (1974)  cheiracterize  these  regions 
as  polar  rain,  polar  showers,  and  polar  squalls.*^  The  polar  showers  occur  during  magnetically 
quiet  conditions  whereas  the  polar  squalls  occur  during  magnetically  disturbed  conditions  and 
are  associated  with  discrete  aurora  in  the  night  sector.  Newell  and  Meng  (1992)  show  that  the 
poleU’  rain  consists  of  low  energy  electrons  of  a  few  hundred  eV  associated  with  the  low  latitude 
boundary  layer  (LLBL),  whereas  the  boimdary  plasma  sheet  (BPS)  is  associated  with  discrete 
aurora  and  the  central  plasma  layer  (CPL)  is  associated  with  difiiise  aurora  on  the  nightside.*® 
Thus,  in  the  day  sector  the  precipitation  has  multiple  sources  depending  on  the  energy  spectrum 
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Figure  4. 
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Figure  S.  Comparison  of  Starkov  Oval  for  Q=9  With  Particle  Precipitation  Ovals  for  Kp=:2  and  4 


of  the  predpitating  particles.  For  the  auroral  clutter  one  is  more  interested  in  the  kmospheric 
irregularities  associated  with  the  aiunral  E-layer  which  is  the  continuous/diflfiise  auroral  oval 
with  a  threshold  of  »>!  keV. 

In  Figure  6  we  reproduce  Figure  (12-12)  from  the  Handbook  of  Geophysics  and  the  Space 

Environment*^  to  emphasize  the  difference  in  thresholds  of  the  respective  instruments.  Figure  6 

shows  the  magnetic  midnight  dependence  of  the  equatorward  boundary  with  Kp  for  electrons 

(DMSP),  auroral  oval  ASCA  photographs,  auroral  oval  from  DMSP  optical  images,  and  the  6300  A 

_  • 

photometer  observations.  The  excellent  match  between  DMSP  electrons  and  OI  6300  A  emission 

shows  that  the  DMSP  electron  sensors  with  lower  thresholds  are  capable  of  measuring  soft 

particle  precipitation  which  produces  6300  A  emission  at  F-layer  altitude.  The  DMSP  optical 

feature:  the  equatorward  boundary  of  the  condnuous/diiBuse  aurora,  (usualfy  observed  aboard 

the  same  satellite  measuring  the  particle  precipitation)  is  coUocated  with  the  particle  precipitation 

boundary  at  very  low  Kp.  and  their  separation  increases  with  the  increase  in  magnetic  activity.  In 

Figure  6.  the  electrons,  the  OI  6300  A  emission,  and  the  DMSP  optical  auroral  boundary  show 

very  good  but  different  dependencies  on  Kp.  Comparatively,  the  ASCA  observations  with  relatively 

poor  resolution  show  more  variation  for  the  Q  dependence  on  I^.  For  ^sterns  operating  at  high 

latitudes,  the  feature  of  Interest  is  the  optical/diffuse/continuous  aurora,  which  originates  at  E- 

layer  altitudes.  The  ionospheric  irregularities  associated  with  the  auroral  E  layer  are  the  source  of 

auroral  clutter,  which  deteriorates  the  radar  performance  at  high  latitude  regions. 

The  least  squared  deviation  fit  to  the  DMSP  diffuse  aurora  on  the  midnight  meridian  in  Figure 
6  yields  the  empirical  equation 


At  =  68.01-1.32 ‘Kp 


(8) 


Dandekar*^  derived  the  empirical  fit  between  the  diffuse  oval  and  Kp 


At=  90-(18.76  +  1.39  •KP+2.68  cosrr*15-18)) 


(9) 


The  corresponding  equation  with  respect  to  Q  is 


At  =  90  -  (20.85  +  1.01  +  2.55  cos(15T  -  9)) 


(10) 


Equations  (8)  and  (9)  are  in  excellent  agreement  with  each  other.  Comparison  of  Eqs.  (6)  and  (9)  shows  that  the  radius 
of  the  particle  oval  is  23.6°  and  the  rate  of  increase  with  Kp  is  1.9°. 
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Kp 

Figure  12- 12.  Ebcamples  of  auroral  equatorward  boundaries  near  midnight. 


Electrons:  Precipitating  electrons  [Gussenhoven  et  al.,  19831® 
ASCA:  All-sl^  camera  photographs  (Feldstein  and  Staiicov, 
1967P®  DMSP;  DMSP  images  ISheehan  and  Carovillano, 
19781*'^  6300  A:  Photometer  {Slater  et  al.,  1980]*®.  Error  bars 
indicate  standard  deviations  for  the  precipitating  electrons  but 
are  also  comparable  for  the  other  measurements. 


Figure  6.  Auroral  Equatorward  Boundaries  Near  Midnight  for  Electrons,  DMSP  Aurora,  ASCA  Aurcna,  and  6300  A 
Photometric  Observations 


whereas  the  corresponding  quantities  for  the  diffuse  auroral  oval  are  21.4°  and  1.4°  respective^. 
Thus,  the  particle  precipitation  oval  Is  2.2°  larger  in  radius  and  e;q>ands  at  a  rate  that  is  0.5°  per 
unit  Kp  larger  than  that  for  the  diffuse  auroral  oval.  The  Starkov  auroral  ovals  being  determined 
from  aU  sl^  camera  ^tems  (ASCA)  are  closer  to  DMSP  optical  features  than  to  particle 
precipitation  features. 

Ekiuations  (6)  and  (7)  assume  that  the  center  of  these  ovals  Is  fixed  and  the  radius  is  a 
frmction  of  the  level  of  magnetic  activity.  The  circles  fit  to  individual  Feldstein-Starkov^  auroral 
ovals  (Holzworth  and  Meng,  1975)  show  that  along  with  the  increase  of  the  radius  of  the  oval,  the 
center  of  the  oval  moves  away  finm  the  CG  pole  with  increasing  magnetic  activity.^^  Thus,  use  of 
Ekjs.  (6)  and  (7)  for  determination  of  the  Kp-Q  relation  will  result  in  inaccurate  Q  values. 

3.  ANALYSIS 

Because  of  the  above  background,  other  possibilities  for  determining  the  Kp-Q  relations  were 
investigated.  The  two  possible  methods  for  such  a  determination  are  1)  a  direct  use  of  the 
observed  Kp  and  Q  indices  and  2)  Kp-  and  Q-dependent  auroral  oval  features  like  various  oval 
boimdaries,  determined  from  the  observed  common  data  sets. 

The  planetary  Kp  index  is  not  available  in  real  time  to  support  the  active  operating  systems. 
Therefore,  the  AWS  generates  its  own  AWS-K  index,  which  is  similar  to  the  Kp  index.  At  present 
the  Q  index  is  available  fipom  a  single  station.  Sodankyla.  This  index  is  not  available  to  AWS  in 
real  time.  Moreover,  Dandekar  (1979)  showed  that  the  Q  index  measured  frnm  a  single  station  is 
not  accurate  at  all  times  because  the  station  cyclical^  moves  closer  and  away  from  the  auroral 
oval  due  to  the  diurnal  rotation  of  the  earth.  The  Q  data  are  reliable  ohty  when  the  station  is 
near  the  auroral  ova  1  around  magnetic  midni^t  Dandekar  (1982)  studied  March  1978-May 
1981  9  and  Kp  and  AWS-K  data.^  for  determining  hourly  dependence  between  Q,  AWS-K  and  Kp. 
The  results  are  presented  in  Table  3,  row  1  for  01-02  UT  when  the  station  Sodanl^a  is  on  the 
magnetic  midnight  meridian.  The  rate  of  increase  of  Q  is  1.4  times  that  of  Kp  and  not  twice  I4>  as 
given  in  Eq.  (3). 

Using  the  larger  data  base  of  the  remaining  22  hours,  Dandekar  used  the  correction  factors*^ 
for  correcting  and  detrending  the  Q  index  data.  The  cumulative  distributions  of  the  corrected 
(detrended)  Q  data  and  the  corresponding  Kp  data  for  the  four  years  are  shown  in  Figures  7 A  and 
7B.  The  figures  show  that  the  distributions  of  both  Kp  and  Q  indices  were  very  similar  in  years 
1969  (close  to  the  maximum  of  the  simspot  c)rcle  20)  and  1971  (approximate  midpoint  of  the 
descending  part  of  sunspot  cycle  20).  The  distributions  in  1965  (closer  to  simspot  minimum  of 
sunspot  cycles  19-20)  and  1974  (towards  the  minimum  of  the  later  part  of  sunspot  cycle  20)  show 
the  largest  differences.  The  data  for  these  years  were  selected  to  represent  a  complete  solar  cycle. 
From  Figures  7A  and  7B,  the  correspionding  values  of  Kp-Q  indices  were  computed  at  every  5 
percent  interval  of  the  population.  These  Kp-Q  values  for  the  four  years  are  shown  in  Figure  8. 
The  least  squared  deviation  straight  line  fit  for  all  the  data  of  each  year  are  shown  in  row  2  of 
Table  3.  On  the  average,  Q  increases  at  the  rate  of  1.56*Kp  and  not  at  twice  Kp  as  given  in  E^.  (3). 
The  table  shows  that  the  results  from  detrended  Q  data  are  in  good  agreement  with  those 
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Figure  7B.  Cumulative  Population  of  Kp  Index  for  Years  1965,  1969,  1971,  and  1974 


restricted  to  the  midnight  meridian  situations  (first  two  blocks  of  Table  3). 

The  preliminary  study  had  shown  that  the  data  should  be  divided  into  two  intervals:  1)  I4>  ^ 
2-i-  and  2)  Kp  >  2+.  The  results  of  the  least  squared  deviation  straight  line  fits  for  these  two  cases 
are  presented  in  rows  3  and  4  of  Table  3  and  in  Figure  8.  When  the  data  are  divided  into  two 
intervals  of  Kp,  the  Kp-Q  dependence  in  these  two  intervals  is  significantly  different  (blocks  3  &  4 
of  Table  3).  One  can  use  either  set  (a  single  range  or  two  Kp  intervals).  The  comparison  of  these 
results  with  Hardy  E^.  (3)  shows  that  for  high  Kp  (>2+)  values  the  slopes  are  in  good  agreement 
but  for  low  Kp  (^+)  values  the  slope  is  about  half  of  that  of  the  Hardy  equation.  The  cutoff  fix>m 
Hardy  matches  with  that  for  low  Kp,  but  for  high  Kp  the  difference  is  quite  significant.  Thus,  in 
general  the  Hardy  equation  does  overestimate  the  Q  values. 

The  equatorward  auroral  oval  boundaries  for  the  midnight  meridian  computed  from  the  Hardy 
equation  are  available  from  the  CRR£^  program."  The  Kp  indices  computed  firom  these 
boundaries  are  shown  with  the  planetary  Kp  indices  for  7687  observations  over  a  six-month 
period  (July-December  1990)  in  Figure  9.  The  middle  curve  shows  the  median  whereas  the  lower 
and  upper  curves  show  the  respective  ±  o  levels.  Up  to  Kp  <  3  both  the  indices  are  in  good 
agreement.  The  Kp  >  3  determined  from  CRRE^  observations  is  slightly  overestimated  with 
increasing  Kp.  Thus,  particle  precipitation  ovals  provide  realty  good  estimates  of  Kp.  It  is  the  Kp-Q 
relation  in  Ek^.  (3)  which  looks  inaccurate. 

The  small  sample  of  AWS-K  data  (in  Figure  2)  is  compared  with  the  Kp  data  in  Figure  10.  The 
figure  shows  that  the  AWS-K  is  overestimated  with  respect  to  Kp  for  Kp  <  3.  The  over-estimation 
of  AWS-K  would  also  result  in  overestimation  of  Q  from  the  Hardy  Ek).  (3).  It  is  essential  to  look  at 
a  larger  AWS-K  data  base  to  see  if  the  over  estimation  of  Kp  in  the  lower  ranges  is  a  persistent 
problem. 

Now  let  us  look  at  the  indirect  determination  of  the  Kp-Q  relation.  Feldstein  and  Starkov 
determined  the  dependence  of  the  three  (poleward,  middle,  and  equatorward)  optical  auroral  oval 
boundaries  from  the  all  slqr  camera  systems  (ASCA)  with  respect  to  Q  and  Kp  indices.  *2  The  fit  to 
the  equatorward  boundary  of  the  Feldstein-Starkov  data  (Figure  2D,  1967)  yields 


X  =  68.93-  1.37  Q 


(11) 


This  equation  is  in  good  agreement  with  that  for  the  diffuse  aurora  observed  aboard  DMSP 
satellites,  shown  in  Figure  6  and  in  E^.  (8)  given  above. 

Feldstein  and  Starkov  graphically  presented  various  auroral  oval  boundaries  as  functions  of 
the  Kp  and  Q  indices.  These  data  are  read  from  the  graphs.  The  detrended  Q  and  Kp 
populations^^  are  used  as  weighting  factors  in  determining  least  squared  deviation  straight  line 
fits  to  the  respective  data  (the  data  points  showing  large  deviations  from  the  straight  line  fit  were 
not  used).  The  respective  relations  tying  latitude  location  to  Kp/Q  were  used  for  determining  the 
Kp-Q  relationships.  These  are  presented  in  Table  4. 

Dandekar  determined  empirical  relations*®  between  the  level  of  magnetic  activity  (Kp  and  Q 
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Table  2.  Parameters  of  the  Starkov  Q  and  Haitfy  Kp  Ovals 
(COMPARISON  OF  STARKOV  (EQ2)  AND  HARDY  (EQ6)  OVALS) 


OVAL 

RADIUS 

OFFSET  FROM 
CGPOLE 

CHANGE  IN  RADIUS 

STARKOV-Q 

(18+5.1)=23.1 

5.1 

.9  PER  UNIT  Q 

HARDY-KP 

(20.27+3.37)=23.6 

3.4 

1.7  PER  UNIT  KP 

Table  3.  Kp-Q  Relation  From  Observed  Kp-Q  Data 


1  EQUATION  FOR  Q*KP 

SOURCE 

REFERENCE 

1  9=1.314K+0.911 

FIGS  2A  &  2B22 

GEOMAG  & 
AERONOMY 

ENG.ED.  7.  1967, 

P49 

Q=1.165K+1.273 

FIGS  3A  &  3C22 

EQUATORWARD  NIGHT 

BOUNDARY 

GEOMAG  & 
AERONOMY 

ENG.ED.  7.  1967. 

P50 

9=1.030K+3.245 

FIGS3A&3C22 

POLEWARD  NIGHT  BOUNDARY 

GEOMAG  & 
AERONOMY 

ENG.ED.  7,  1967, 

P50 

9=1.716K-4.413 

FIGS  4A  &  4C22 

POLEWARD  DAY  BOUNDARY 

GEOMAG  & 
AERONOMY 

ENG.ED.  7.  1967. 

P50 

Q=2.167K-2.058 

FIGS  4A  &  4C22 

EQUATORWARD  DAY  BOUNDARY 

GEOMAG  &  1 

AERONOMY 

ENG.ED.  7.  1967. 

P50 

9=1.506K+0.625 

FIGS  5 

OVAL  WIDTH 

GEOMAG  & 
AERONOMY  ENG.ED. 

7.  1967.  P51 

9=1.376K-1.9 

Q=1.333K-2.0 

N.  HEMISPHERE  DMSP  DATA® 

S.  HEMISPHERE  DMSP  DATA® 

AFGL-TR- 

79-0010 
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KP  DETERMINED  FROM  CRRES  DPlTfi 


OBSERVED  KP  VRLUES 


PLRNETRRY  KP  INDEX 


Figure  9.  Comparison  of  Auroral  Kp  Index  With  the  Planetary  Kp  Index  From  CRRES  Data  of  July-December  1990 
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Figure  10.  Comparison  of  the  AWS  -K  Index  With  the  Planetary  Kp  Index 


Table  4.  Kp-Q  Relation  From  Oval  Parameters 


EQUATION 

REFERENCE 

REMARKS 

9=1.37  'KP-l.? 

(01,02  un 

DANDEKAR( 1982)5 

AFGL-'m-82-0010 

SODANKYLA.FINLAND 

67.4  N.  26.6  E 

9=1.245 'KP -0.5 

9=1.519*  KP-l. 1 

9=1.581  *  KP-1.0 

9=1.823*  KP -2.0 

9=1.558  *KP -1.1 

PRESENT  STUDY 

1965 

1969 

1971 

1974 

ALL  FOUR  YEARS 

9=0.962  *  KP  -0.2 

9=0.989  *  KP  -0.4 

9=1.007  *  KP -0.3 

9=1.065  *  KP -0.6 

9=0.964  *  KP  -0.3 

PRESENT  STUDY 

1965 

1969 

1971 

1974 

ALL  FOUR  YEARS 

FOR  KP  <  2+ 

9=1.871  *  KP-2.4 

9=1.950  *  KP -2.4 

9=2.000  *  KP  -2.4 

9=2. 189  *KP -3.4 

9=2.040  *  KP  -2.7 

PRESENT  STUDY 

1965 

1969 

1971 

1974 

ALL  FOUR  YEARS 

FOR  KP  >  2+ 

9=1.125*  KP -0.7 

9=1.497*  KP-1.0 

9=1.302  *  KP -0.7 

9=1.545  *KP -1.4 

9=1.367*  KP -0.9 

DANDEKAR( 1982)5 
AFGI^TR-82-0010 
(ADAl  18734) 

1965 

1969  1 

1971  1 

1974  1 

ALL  FOUR  YEARS  1 
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(not  detrended))  and  the  equatorward  boundary  of  the  optical  continuous/difiuse  of  both 
northern  and  southern  hemispheric  aurora  observed  aboard  the  DMSP  satellites.  The  results  are 
shown  in  Table  4. 

Thus,  in  general,  results  in  Table  4  show  that  Q  increases  approximately  at  the  rate  of  1.5*Kp 
and  not  at  twice  the  rate  of  Kp  as  determined  by  Hardy  in  Eki-  (3). 

4.  CONCLUSIONS 

The  particle  precipitation  auroral  ovals  do  provide  a  good  real  time  estimate  of  the  Kp  index. 
In  the  above  we  see  that  matching  the  particle  precipitation  auroral  ovals  with  the  Starkov  optical 
auroral  OArals  leads  to  an  inaccurate  determination  of  the  Kp-Q  relationship.  It  would  be  nice  to 
have  a  midnight  meridian  Q  index  to  compare  with  the  planetary  Kp  index  for  determining  the 
Kp-Q  relation.  However  a  reasonably  direct  relationship  between  Kp  and  Q  is  obtained  by 
correcting  the  observed  single  station  Q  data  for  the  diurnal  rotation  of  the  earth  under  the 
auroral  oval.  For  a  better  accuracy  in  the  determination  of  the  Q  index,  the  Kp  range  is  divided  in 
two  intervals.  The  improved,  more  precise  Q  values  obtained  from  this  procedure  will  improve  the 
prediction  of  the  auroral  oval  phenomena  that  affect  the  performance  of  the  real  time  operational 
systems  in  high  latitude  regions. 
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